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Abstract 

We demonstrate the long range strong coupling of magnetostatic modes in spatially separated 
ferromagnets mediated by a microwave frequency cavity. Two spheres of yttrium iron garnet are 
embedded in the cavity and their magnetostatic modes probed using a dispersive measurement 
technique. We find they are strongly coupled to each other even when detuned from the cavity 
modes. We investigate the dependence of the magnet-magnet coupling on the cavity detuning A, 
and find a 1/A dependence also characteristic of cavity-coupled superconducting qubits. Dark 
states of the coupled magnetostatic modes of the system are observed, and ascribed to mismatches 
between the symmetries of the modes and the drive field. 
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There are several mechanisms by which spin angular momenta can couple, producing new 
total spin eigenstates. In most circumstances, the dominant coupling is either the dipole 
interaction, where the magnetic dipole moments interact through the electromagnetic field, 
or exchange coupling through a combination of the wavefunction symmetry and electrostatic 
interactions. The coupling of spins does not always occur directly; in many cases indirect 
mechanisms play an important role. The additional element in an exchange-type interaction 
can be a localised electron as in super-exchange, or an itinerant carrier as in the RKKY 


interaction[l|, which enables magnetic coupling between two ferromagnets through a thin 
electrically conducting paramagnetic spacer layerj2-4|. The dipole interaction can also be 
involved in indirect coupling, such as in the J-coupling of nuclear spins mediated by a 
combination of the local dipole coupling between the nuclear spin and electrons and exchange 
coupling of the electrons on the separate nuclei P 


An alternative indirect coup 
context of cavity and circuit 


)hng between spin-like objects has been demonstrated in the 
3, ISl] QED; in this case the additional element is a low-loss 
resonator. It is not necessary for the psuedospins to be resonant with the cavity. Instead 
the spins and cavity can be significantly detuned from one another, in what is known as the 
dispersive regime. Here, an interaction between the pseudospins and the cavity via the local 
dipole coupling leads to an effective coupling between the individual pseudospins mediated 
by virtual photons[7|, l9|. This approach can be used to couple together quantum systems 
such as qubits in a controlled way at distances far beyond that achieved by the qubit’s dipole 


alone 
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12| . The versatility of this approach is demonstrated by the varietj^f a 
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ternative 


16| . double 


systems which can be coupled to the resonator, including spin ensembles 
quantum dotsjn, 18| and hybrid systems 3, 

Here we apply a dispersive technique to demonstrate long-range coupling of the macro¬ 
scopic magnetic moments of two ferrimagnetic yttrium iron garnet (YIG) spheres mediated 
by an electromagnetic cavity. Due to its low damping YIG is an important material for mi¬ 
crowave coimonents such as tunable filters and couplers 211] , as well as current research into 
spintronics 22N24| . Stroiig coupling between YIG magnetostatic modes and microwave cavi¬ 
ties is readily attainable 25l-l28|. and has been exploited in recent work 29| towards coherent 
control of single magnon states. In such systems changes in longitudinal magnetisation shift 
the cavity mode frequency j^. yielding a dispersive measurement of ferromagnetic resonance 
(FMR). 
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In our experiments the two 1 mm diameter YIG spheres are placed within the dielectric of 
a co-axial transmission line cavity, made by breaking the inner of a 3.5 mm semirigid coaxial 
cable in two places. They are positioned at the antinodes of the second harmonic resonance 
(a;2/27r = 7.15 GHz, loss rate k/2tt = 45 MHz) (Fig. la). A magnetic held is applied 
in the coaxial direction, consisting of a uniform held, Hq, and a diherential held, ±6H/2, 
local to each sphere. Hq is sufficiently large to saturate the magnetisation of the YIG. In 
addition to the uniform magnetostatic mode of linewidth z//27r 5 MHz, a single domain 

ferromagnetic sphere hosts a spectrum of non-uniform magnetostatic modes which lie in a 
frequency band around the uniform modejM, 321. We focus our attention on the uniform 
mode, corresponding to the magnetisation processing in phase throughout the sphere; this 
couples most strongly to the resonator held. 

In order to characterise the strength of the resonant magnet-cavity coupling we measure 
the transmission of a -10 dBm probe tone of frequency LJp/27r through the cavity. Initially 
we detune the modes of the two spheres from each other by ~ 450 MHz {SH ^ 16 mT) 
and measure the cavity transmission as a function of Up and Hq around its second harmonic 
(Fig. lb). The FMR frequencies of the spectrum of magnetostatic modes are seen to come 
separately into resonance with the microwave cavity and the avoided crossin gs s how that 


the uniform modes of the YIG spheres are strongly coupled to the cavity 


25 
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29| . For this 

cavity mode we determine a magnet-cavity coupling of gi^^/2TT k, 150 MHz for both spheres. 


and we hnd the coupling to higher order magnetostatic modes to be weaker 28(] . Here, 
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100, conhrming that the system is in the strong coupling regime. Similar behaviour 
(not shown) occurs around the fundamental mode of the cavity, u:i/2tt = 3.55 GHz, with 
gu:i/2TT K. 80 MHz. 

Having characterised the resonant coupling between the magnets and the cavity, we now 
move into the dispersive regime by adjusting Hq such that the magnetostatic modes are 
signihcantly detuned from the cavity modes. In order to probe the magnetisation dynamics 


of the magnets we use a dispersive measurement technique 


33 


34j . in which we measure the 


phase 0 of the transmitted probe signal at a resonant frequency of the cavity, Up = 
while applying a second drive tone ujd of power 27 dBm to excite FMR. The measured 
change in phase when FMR is driven is due to the reduction in the sum of the longitudinal 
components of magnetisation of the spheres 30j. Fig. lb (monochrome data) shows the 
dispersive measurement of the uniform FMR modes of each individual magnetic elements 
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FIG. 1. Experimental scheme, (a) Experimental setup. Two YIG spheres are positioned at the 
magnetic field antinodes of the second harmonic of a transmission line cavity. Two sources, at 
ojp and ujd, are coupled into the cavity. The transmitted amplitude and dispersive phase shift at 
iVp is measured by homodyne detection. A global field, Hq, is applied to align the magnetisation 
of the spheres in the propagation direction of the cavity, and tune the ferromagnetic resonance 
to be off-resonance with the cavity modes. The field at each sphere is adjusted by an amount 
±5Hj2 using a local coil wrapped around the cavity, (b) Response of the cavity as measured by 
transmission measurements of |S21| (with ujp between 6.625 GHz and 7.5 GHz and the source at 
(jJd turned off - coloured data) and dispersive measurements (with Up = uji, and ujd between 5.5 
GHz and 6.5 GHz, and 7.625 GHz and 9 GHz - greyscale data). 


above and below 002 , with them still detuned from each other by SH = 16 mT. Using this 
technique, we are able to measure FMR far detuned from the cavity modes. 

In order to investigate the magnet-magnet coupling, Hq is hxed such that the FMR 
frequencies of the uncoupled magnets (wpi, 00 ^ 2 ) are ~ 0.8 GHz below a; 2 , and Up is set to Ui 
to avoid resonant interactions between higher order FMR modes and the measurement cavity 
mode. We now attempt to bring the two magnetic modes through resonance by varying 6H 
from —6 mT to -1-6 mT. An avoided crossing is observed (Fig. 2a) demonstrating coupling 
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FIG. 2. Magnet-magnet coupling with and without a cavity, at Hq = 200 mT. (a) Dispersively 
measured spectroscopy of the magnetostatic modes in the presence of a cavity close to and below 
0 J 2 . An anticrossing between the modes is observed. The additional structure is due to weaker 
coupling to non-uniform magnetostatic modes, (b) Transmission amplitude |S21| of the system 
with the coupling gaps closed, and therefore no cavity present. No anticrossing is seen, with the 
magnetostatic linewidth of 5 MHz placing an upper bound on the magnet-magnet coupling. 

of the magnetisation dynamics of the two spheres. The magnitude of the coupling is given 
by the frequency splitting of the modes at 6H = 0, giving a value of 2J/2tt = 87 MHz for 
this cavity detuning. 

We hrst exclude the possibility that this frequency splitting is simply due to the free 
space dipole coupling between the magnets. By closing the coupling gaps at both ends of 
the cavity, and reverting it back to a simple coaxial cable, the cavity modes are eliminated 
but the spatial separation of the magnets, and therefore the direct dipolar coupling, remains 
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the same. In Fig. 2b we show |S21| of the resulting transmission line as a function of 6H, and 
observe the magnetostatic modes of the two spheres come into resonance with each other. 
This demonstrates that the splitting requires the cavity to be present, and is therefore not 
due to direct dipolar coupling. In addition, a simple calculation of the dipolar coupling 
can be made by linearizing the Landau-Lifshitz equation for the two individual magnets 
and including the dipolar stray held from each at the other. The harmonic solutions of the 


resulting coupled equations give a coupling rate of 2 J = 7 


27 r|r| 




For YIG with a saturation 


magnetization of rus = 140 kA m“^ and a magnet-magnet separation of r = 14 mm, and 
setting the gyromagnetic ratio to 7/27r = 28 GHz T“^, we obtain a value of ~ 200 kHz, 
much less than the coupling we observe. 

The strong magnet-magnet coupling we observe in the presence of the cavity is analogous 
to that between qubits embedded in a microwave cavity lOj , which is generally described 
using the Jaynes-Gummings Hamiltonian 35(]. Here, we adopt a quantum mechanical model 
of the Kittel mode of each magnet in the macrospin approximation 36l-l38| . The magnetisa¬ 
tion of the hrst (second) magnet is modelled by a large quantum angular momentum 51 ( 2 ). 
The full Hamiltonian for a system with two magnetic spheres in the cavity is then 


n = j{Ho + 6H/2)Su + 7(^0 - SH/2)S2, 

+'y6Hy{a + a^)(5 lx T S2x) + TiojcCi^a. 

This is the sum of the Zeeman energy of the two magnets in the total magnetic held 
including the cavity mode with r.m.s vacuum magnetic held 5H, and the photon energy of 
the cavity held with frequency ojc and lowering operator a. The YIG spheres are located at 
antinodes of the cavity held and the a.c. magnetic held lies along the x-direction, with the 
sign in the third term depending on the relative phase of a.c. held at the two spheres. 

We now move into an interaction picture with respect to the FMR resonance frequency 
of the Kittel mode. For simplicity we consider 5H = 0, and we dehne A = Uc — qhfo to be 
the detuning of the cavity resonance from the FMR frequency. We also make the rotating 
wave approximation and obtain the interaction picture Hamiltonian 


Hint h^o[®(5l+ i ^2+ ) + at(5i_ ± 52-)] + hAata, (2) 

where we have dehned the single spin coupling frequency go = 'jSHy. 
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In the dispersive limit where A 3> one can readily obtain, by applying second order 
degenerate perturbation theory to "Hint, the effective Hamiltonian 


2 ( 7 ^ 

•Heff = ± 

+ + ^(*S'l2 + 5'22)] 

0^2 

+ J^a)a{Si, + 52,) + Aa^a. 


( 3 ) 


This is perfectly analogous to the more familiar case of two qubits in a single cavity j^. In 
the first term we find a coupling of the transverse components of the two magnetic moments, 
as observed in the experiment. The second term describes both linear and nonlinear shifts 
in the FMR resonance frequency and the third term is the dispersive coupling of the Kittel 
modes to the cavity which permits our dispersive measurement. The final term is the cavity 
detuning. 

In the case where each YIG sphere is highly polarised it is usual to analyse this model 
in terms of the Holstein-Primakoff transformation j^|, in which the magnetisation of the 
YIG sphere can be described in terms of a harmonic oscillator mode with lowering operator 
b for which 5, = h(Yeff/2 — b^b), where is the number of unpaired spins per magnet. 
Writing the first term of equation ([3]) in terms of b we find the coupling Hamiltonian Tic = 
±^{SixS 2 x + SiyS 2 y) ~ hJib\bi + b\b2)/2 where J = ±2Yefj5fg/A. The expected normal 
mode splitting resulting from the coupling J is dependent on the total coupling = g^y/N^ 
of the magnetostatic mode to the nth cavity mode and inversely proportional to the detuning 
A. 

In order to verify that this is the correct coupling mechanism, we now study the depen¬ 
dence of J on A. We measure J over a large range of FMR frequencies on both sides of the 
072 mode of the cavity (Fig. 3), whilst remaining in the dispersive regime. Although in Eqn 
[3] we have explicitly discussed only the effect of a single cavity mode it is straightforward to 
add other cavity modes to the model to obtain the formula for the observed normal mode 
splitting. In the measurement range there are three modes that play an important role in 
the coupling. When the coupling is dominated by a single mode, we observe the predicted 
1/A dependence, as previously observed in circuit QED[40|. We fit to this data the sum of 
the couplings due to the three modes, with corresponding magnet-cavity mode couplings, 
g^^, as free parameters. Our model is in good agreement with the data, and we extract 
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FIG. 3. Absolute magnet-magnet coupling, 2| J|/27r, as a function of FMR frequency of the uncou¬ 
pled modes; this is determined by Hq. Blue squares are measured with Wp = cui, and red circles 


are measured with ujp = U 2 - Orange regions denote the resonant strong coupling regime 


33l |. where 


A < The curve is from the model described in the text summed over the lowest three cav¬ 


ity modes, with values of magnet-cavity couplings as fitting parameters. We extract couplings of 
g^^j2TT = 88 MHz, = 177 MHz and = 83 MHz, in approximate agreement with the 

coupling deduced from transmission measurements. 


values of = 88 MHz, g^^/2Ti = 177 MHz and g^^l2T\: = 83 MHz, consistent with the 

resonant couplings obtained from transmission measurements. 


The sign of J is governed by the sign of the detuning and the spatial symmetry of the 
coupling modes, and determines whether the lowest energy state of the coupled system is 
symmetric {Si^ = S 2 x) or antisymmetric {Six = —S 2 x)- These modes correspond to in 
phase and out of phase precession of the two magnetisations. In Figs 4a and 4b, we show 
two measured anticrossings above and below UJ 2 - In both cases one of the eigenmodes cannot 


be driven; these dark states are similar to those seen in coupled superconducting qubits 


4Q|. 


They are due to the microwave drive in the cavity forming a standing wave; at the sites 
of the two magnets the amplitude is in general different, and the phase difference can be 
either 0 or tt. If this does not match the the symmetries of the coupled FMR mode, then 
the mode cannot be driven and is dark. In Figs 4c and 4d we show the predicted response of 
the system to a drive tone, based on the calculated form of the drive held in the cavity, the 
symmetry of the coupled mode as a function of 6H and the measured dispersive lineshape 
of the uncoupled modes. The visibility of the modes in the calculated response agrees well 
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FIG. 4. Dark states and eigenmode symmetries, (a) Dispersively measured spectroscopy of the 
anticrossing between magnon modes below UJ 2 and (b) above 0 J 2 - (c) and (d) Modelled spectroscopy 
for a and b, taking into account detuning from the cavity modes and the frequency dependent drive 
symmetry, but not including non-uniform magnetostatic modes. 


with our data. 


In conclusion, by measuring avoided crossings of magnetostatic modes we have shown 
that spatially separated magnetic moments may be passively coupled over a long range 
via dispersive coupling to the modes of an electromagnetic cavity. The coupling can be 
understood within the framework of circuit QED, putting magnets on a similar basis to 
qubits and atoms in cavities. This suggests a route towards coupling magnets to a variety 
of physical objects including mechanical oscillators and superconducting qubits using the 
dispersive strong coupling regime of QED, a different approach to that demonstrated by 
Tabuchi et al MJ- Such an approach might be used in coherent magnetic metamaterials, 
or to phase-lock many spatially separated magnetic oscillators, such as those in spin-torque 


nano-oscillators 


42| . In such a scheme, the dispersive nature of the interaction means the 


linewidth of the coupled oscillators would not be limited by the quality of the cavity. Indeed, 
the linewidth of the cavity can be greater than the coupling rate between the magnets 10 1. 
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